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Abstract 

The LEPS/SPring-8 experiment made a comprehensive measurement of the spin-density matrix 
elements for ^yp — > <fip, 'yd — > (fipn and — > (fid at forward production angles. A linearly polarized 
photon beam at E 1 = 1.6 — 2.4 GeV was used for the production of cfi mesons. The natural-parity 
Pomeron exchange processes remain dominant near threshold. The unnatural-parity processes of 
pseudoscalar exchange is visible in the production from the nucleons but is greatly reduced in 
the coherent production from deuterons. There is no strong E 1 dependence, but there is some 
dependence on momentum transfer. A small but finite value of the spin-density matrix elements 
reflecting helicity-nonconserving amplitudes in the t-channel is observed. 
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I. INTRODUCTION 



The study of photoproduction of light vector mesons (p, u, and 0) plays an important 
role in hadron physics to understand the nonperturbative aspect of QCD. At high energies 
(W > 10 GeV), Pomeron exchange in the t channel describes well the photoproduction 
mechanism in the framework of vector-meson dominance. The Pomeron is introduced in 
the Regge theory for describing high-energy hadron scattering, and it is generally believed 
to originate from multi gluon exchange processes. At low energies near threshold (W ~ 2 
GeV), the reactions are open to the exchange of mesons and baryons, where possible "missing 
resonances" could be involved. In particular, the gluon dynamics dominates in the 0-meson 
production process because the quark exchange is suppressed by the Okubo-Zweig-Iizuka 
(OZI) rule owing to the major ss content of the meson. Therefore, the photoproduction of 
mesons at low energies provides a unique opportunity to explore the behavior of Pomeron 
exchange and to study exotic hadronic interactions mediated by multi gluon exchanges. 
In fact, 0-meson photoproduction from protons and deuterons has been measured in the 
Laser Electron Photon beamline at SPring-8 (LEPS) [l-4], the CEBAF Large Acceptance 
Spectrometer (CLAS) at Jefferson Lab 5j47|], and the Spectrometer Arrangement for Photon 
induced Reactions (SAPHIR) at Bonn 8]. 

Spin observables are known to be a powerful tool to obtain further insights into the rel- 
evant reaction mechanisms. With the use of linearly polarized photons, the decay angular 
distribution of the vector meson can be expressed in terms of the real or imaginary part of 
nine spin-density matrix elements (pf ) |9[]. These matrix elements are sensitive to the under- 



lying reaction mechanism 



10j-[l4|. For example, the matrix element p\_ l reflects the asym- 



metry in the contribution from natural-parity and unnatural-parity exchange [12||. There 



were some measurements of decay asymmetry [15j and spin-density matrix elements 
of production from free protons at high energies (2.8, 4.7, 9.3 and 20—40 GeV). The results 
were consistent with s-channel helicity conservation. Judging from the fact that p\_ l 0.5, 
one can conclude that the meson is produced predominately by natural-parity Pomeron 
exchange at high energies. 

At E 1 = 1.6 — 2.3 GeV, the photoproduction of mesons from free protons in the forward 
direction [l| can be mostly described by the Pomeron and (n, rf) exchange in the t channel. 
A value of 0.2 was found for p\_ x , which deviates from 0.5, the limit corresponding to 



pure natural-parity exchange. It suggests a nonnegligible contribution of unnatural-parity 
exchange processes, in contrast to complete dominance of the Pomeron exchange at high 
energies. 

In the measurement of coherent <fi production from deuterons [2|, it is reported that p\_ 1 
clearly becomes close to 0.5. This suggests that the dominant unnatural-parity component, 
isovector tt exchange, is forbidden in the coupling to the isoscalar deuteron target (3, [jjj . 

The nuclear transparency ratio of <fi photoproduction for deuterons shows a large sup- 
pression of incoherent production 3] and is consistent with the A dependence of the ratio 
observed for nuclear targets ^]. Also, p\_ x was observed to be slightly larger in incoher- 
ent production from deuterons, compared with production from free protons. It coincides 
with destructive interference between isovector tt and isoscalar rj exchange amplitudes in the 
7n — > <fin reaction [10, 111]. However, there is no observation of a large isospin asymmetry 
of (fi production from nucleons. The decrease in the 0-meson yields, scaled as number of nu- 
cleons in the production from deuterons, cannot be adequately explained in terms of isospin 
asymmetry. 

A narrow bump structure was found around E 1 = 2 GeV in 0-meson production cross 
sections from nucleons [l] and the origin of this structure is not yet understood. Naively, 
it is speculated to result from the appearance of (tt,t]) exchange near threshold. However, 
the measured decay angular distributions, which in principle reflect the relative weights of 
natural-parity and unnatural-parity processes, do not vary appreciably across the bump re- 
gion. This structure cannot be described by a conventional model with the Pomeron and 
pseudoscalar exchange, where a monotonic energy dependence is predicted Various the- 
oretical interpretations of this structure are proposed: the interference between unnatural- 
parity tt and 77 exchanges [l^] , the interference of the isovector scalar ao meson with Pomeron 



exchange the coupled-channel effect of JCA(1520) and <fiN 



and the existence of N* 



resonances with a large ss component 



19] . Nonetheless, a satisfactory description of the 



bump structure of 0-meson photoproduction is still lacking. Comprehensive information on 
spin-density matrix elements will help clarify the situation. 

Previously, only limited numbers of spin-density matrix elements were determined from 
one-dimensional decay angular distributions [l-3, 6-8]. This paper presents the results of 
a complete set of spin-density matrix elements in the 779 — > <fip, 'yd — > (ppn, and 7c? — >• 
<fid reactions measured in the LEPS/SPring-8 experiment. The angular distributions were 



measured via the charged-kaon decay mode of the meson. In Section |TT] we describe the 
formulation of spin-density matrix elements and the extended maximum likelihood fit used 
for the determination. The experimental setup and the analysis details are introduced in 
Sec. IIII1 Sec. IIVI shows the spin-density matrix elements of photoproduction of <p mesons 
from protons and deuterons. The interpretation of the data and discussion are given in 
Sec. [V] Finally, Sec. ED provides a summary. 



II. SPIN-DENSITY MATRIX ELEMENTS AND EXTENDED MAXIMUM LIKE- 
LIHOOD 

A. Decay angular distribution 

The angular distribution, W, of K + decaying from mesons produced with a linearly 
polarized photon beam can be expressed as follows jjj: 

W(n,$;p) =W (cos 9, <p,$;p) 
=W° (cos 9, (p;p) 

- P 7 cos2$W rl (cosfl,<p;p) 
-P 7 sin2W 2 (cos#,y?;p) 



(1) 



^ = ^[^(l-^ ) + ^(3pgo-l)cos^ 

-V2Rep% sin 2$ cosy? — p ( l_ 1 sin 2 #cos2<p], (2) 

W 1 = — [p} x sin 2 9 + pL cos 2 9 
Aix 

sin26>cos<p — p}_ : sin 2 9cos2ip\, (3) 
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W 2 = — [v^Imp^ sin 29 sin <p 

4:71 

+Imp 2 _ 1 sin 2 6'sin2<p] (4) 

Here, p represents the measurable parts of the nine independent spin-density matrix elements 
pfj, 9 and <p are the polar and azimuthal angles of decay particles in the rest frame of vector 
mesons, $ is the azimuthal angle of the photon electric polarization vector with respect to 
the production plane of vector mesons, and P 1 is the degree of linear polarization of incident 
photons. 



Conventionally, there are three different choices of the quantization axis z' for the decay 
reference system J9]: the helicity system with z' opposite to the velocity of the recoiling 
nucleons in the vector-meson rest frame, the Gottfried- Jackson system with z' parallel to 
the momentum of the photon in the vector- meson rest frame, and the Adair system with z' 
parallel to the photon momentum in the overall center-of-mass (CM) system. 

The physical property of the production mechanism is simplest when illustrated in: 
the helicity system for s-channel helicity conservation, the Gottfried- Jackson system for 
^-channel helicity conservation with no absorption, and the Adair system for the spin in- 
dependence in the overall CM system {20]. In addition, the Gottfried- Jackson system has 
the advantage that some spin-density matrix elements work as a measure of the asymm etry 
between processes with natural-parity and unnatural-parity exchanges in the t channel 11 ]. 
In general, the values of spin-density matrix elements depend on the system chosen. Never- 
theless, the difference among the three systems becomes small at very forward angles. The 
results in the following are presented in all three systems for easy comparison. 

The contribution to the cross section from natural-parity and unnatural-parity exchanges 
in the t channel, a N and a u , can be determined from the density matrix elements. In the 
case of helicity-conserving exchanges, all spin-density matrix elements become zero, except 
for p\_ x and Imp^_ x [12|. That p\_ 1 = — \mp\_ x = +0.5 (—0.5) corresponds to the case with 
pure natural-parity (unnatural-parity) exchange. The matrix element p\_ 1 directly relates to 
the asymmetry in the contributions from natural-parity and unnatural-parity exchanges 12 1 
and is expressed as 



1 l ° N - aU <<\ 



In principle, nonzero values for spin-density matrix elements other than p\_ l and lm.p\_ 1 



indicate the nonconventional OZI-evadine 



)rocesses such as the ss knockout 



21] or non- 



diffractive baryon resonance production 12j, |22j. Those processes are commonly expected 
to become important at large transferred-momentum (|i|) regions for photoproduction. 

The matrix element p(j reflects the strength of the single-spin-flip exchange amplitude 
resulting from the components other than the Pomeron exchange. For evaluation of the 
relative contribution of the unnatural-parity exchange at small one need take into account 
Poo 3 as follows: 

1 a-2 P i 1 -p° 00 ). (6) 



a N _j_ a U _|_ a spin-flip 2 



A finite value of the spin-density matrix element p^-i comes from the amplitude responsi- 
ble for the double-spin transitions where the helicity of the (f) meson differs from the photon 
helicity by two units: A 7 = ±1 — > A^ = =Fl. In many theoretical models for 0-meson pho- 
toproduction, such as scalar, pseudoscalar t-channel exchange, and the original Donnachie- 
Landshoff Pomeron exchange model based on the Pomeron-isoscalar photon identity, these 
transitions are forbidden, and p°_ 1 is exactly zero. In the modified Donnachie-Landshoff 
model, p5-i m the Gottfried- Jackson system could be nonzero owing to the spin-orbital 



interaction inherent to the two-gluon exchange in the t channel ll|, |23 |. 



B. Extended maximum likelihood fit 

Experimentally we measure the angular distribution W(Q, $; p) for charged kaons de- 
caying from photoproduction events of the meson. A few spin-density matrix elements 
can be extracted from the one- dimensional angular distributions of cos 9, cos ip, cos and 



cos(</9 — $) 12|. In this work, a binless maximum likelihood fit is applied to determine the 
nine spin-density matrix elements simultaneously. 

A likelihood function with a perfect detection efficiency, C, is defined as: 



N 



i=1 (7) 
fj / $; p)dVld<$> 

with Pi the normalized likelihood for each event and W(Q, <&; p) the probability density 
function. The logarithmic likelihood function — In £ is minimized by the CERNLIB MINUIT 
package Q. 

If the detector efficiency is not perfect, the probability density function is modified as 

= WiSk^p) = WjQi^up) 

1 fw{n,$;p)r){n)dnd<$> w(p) ' 

By taking into account the efficiency rj(Q), the evaluation of the normalization factor W 
with a variation of p in the minimizing process, is technically challenging. This issue can 
be solved by representing the angular distribution W(Q, $; p) with a set of orthogonal bases 
O lm [25|: 

W(n,*-p) = yj^a /m (p)O im (n,$), (9) 

Lm 



TABLE I. Angular coefficients a/ m and their corresponding bases O 



Lm a lm O lm 

O' 1 vW§ Y o°(^ 



0,2 ^^(3^-1) Y 2*(^k) 

0,3 -yfayfl Re(p? ) ^(-y 2 l- + F 2 - 1 *)(^=) 

0,4 -V^y^pt! ^(r 2 2 * + *7 2 *)(^) 

1.1 -VZP-yJkUo + 2Pri) ^o°(^ cos(2$)) 

1.2 -0FP 7 yj (/9 i o F 2 °*(-^cos(2$)) 

1.3 V^P^ 5 p\o ^(- Y 2* + Y 2 U )(^ cos(2$)) 

1.4 v^V^i-i + >T 2 *)(^ cos(2d>)) 

2.1 -^5FP 7) /f Im(pf ) ^(tf* + y 2 - U )(^ « in ( 2$ )) 

2.2 -0FP 7V /flm(P?-i) i^(^-y 2 " 2 *)(^?Bin(2$)) 

where the orthogonal conditions for each base are satisfied: 

J *lm l'm> dQ = 5ll ,5 mm/ , (10) 

and the coefficient for each basis ai m is 

J 0* lm W(n, $; p) rffi d$ = \j^a lm {p) (11) 

The explicit representations of 0' m in terms of spherical harmonics, and the coefficients 
a; m 's are given in Table [fl 

The angular distribution of the detector efficiency T](p.) can also be expanded by the same 
orthogonal basis O lm : 

r ] (Q) = Y,bimO lm (Q). (12) 
In a kinematic bin, the angular moment of experimental acceptance, bi m , is evaluated as 



blm = \l~kJ 0/m( ^' ^(tydnd® (13) 

^ accepted 

v E ^ m (a,$,), (14) 

* generated 



i=l 

where the numerator and denominator sum over the accepted and isotropically generated 
Monte Carlo (MC) events, respectively. 
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The evaluation of the normalization factor W(p) is simplified as follows: 

w(p) = Jw(n,$]p)ri(n)dnd$ (15) 

= ^2 ai m (p)b lm . (16) 

l,m 

Furthermore, the restriction that the probability distribution P is normalized to 1 could be 
relaxed in an "extended maximum likelihood fit" 26|. The integral of the unnormalized 
probability density function V represents the total number of predicted events, N, under 
the assumption of a Poisson variation for the measured number of events, N. 
The extended likelihood function can be written as 



1=1 



V = yW{n " $ " B* (18) 
1 fyW(tt,<£;p)ri(tt)dttd® 1 ' 

N = J yW{tt,$;p)ri(tt)dttd$ = yW(p), (19) 

where y is the actual yield. By neglecting the terms that do not depend on the fit parameters, 
the likelihood function to be minimized can be rewritten as follows: 

N 

- In C(p, N) = -J2 Hyi(n i; p)\ + N (20) 

i=i 

N 

= -J2 ln[iV/(O i; p)/W(p))] + N. (21) 

i=l 

Now besides p, the expected number of events, N, becomes an additional parameter in the 
minimization, and it should turn out to be the number of events to be fit, N, in the final 
result. This fact could be utilized to validate the success of the minimization procedure. 



III. EXPERIMENT AND ANALYSIS 



The experiment was carried out at the LEPS/SPring-8 facility using a linearly polarized 
photon beam produced by backward Compton scattering of the Ar laser from 8-GeV elec- 
trons in the storage ring of SPring-8. Photons in the energy range of 1.5—2.4 GeV were 
tagged by detecting recoil electrons. The photon beam with an intensity of ~10 6 /s was 
directed onto liquid hydrogen or deuterium targets inside a 15-cm-long target cell. The 
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direction of linear polarization was controlled vertically or horizontally by using a half-wave 
plate for the laser with a polarization of nearly 100%. Charged particles emitted from the 
interaction points of photons with a target were detected at forward angles in the LEPS spec- 
trometer. The spectrometer consisted of a start counter, a silica-aerogel Cerenkov counter, 
a silicon vertex detector, a dipole magnet, three multiwire drift chambers, and a time-of- 
flight wall. The angular coverage of the spectrometer was about 0.4 rad and 0.2 rad in the 
horizontal and vertical directions, respectively. Particle identification was made by mass re- 
construction using the measured time of flight and momentum. For more details concerning 
the detector configuration and the quality of particle identification, see Ref. {23] . In the 
present work, the integrated number of tagged photons reached at 5.6 xlO 12 (4.6xl0 12 ) for 
the hydrogen (deuterium) runs. 

Events with both K + and K~ tracks detected were selected. The spectra of a Dalitz 
plot of the K + K~p final state, K + K~ invariant mass, missing mass and K~p invariant 
mass for those events from hydrogen and deuterium targets are shown Fig. [TJ and Fig. 121 
respectively. A clear peak in the nominal mass of meson, 1.019 GeV, is seen in Fig. UJb) 
and Fig. [2](b). The missing mass spectrum of the events from hydrogen shows a peak 
around the proton mass as seen in Fig. [T](c) while more complex structure is observed in the 
spectrum of the production from deuterium (see Fig. [2]Jc)). In the missing mass spectrum, 
assuming the whole deuteron as the target (MM ( j(7, 0)), events of coherent production, 
7<i — > (pd, peak at the deuteron mass of 1.875 GeV/c 2 whereas incoherent events, 7<i — > (ppn, 
are distributed at relatively higher mass. This missing-mass spectrum is nicely reproduced 
by MC simulations of coherent and incoherent production processes. The MC simulation 
takes into account experimental parameters such as geometrical acceptance, energy and 
momentum resolutions, and the efficiency of detectors. The effects of Fermi motion, along 
with off-shell aspects of target nucleons inside deuterium and final-state interaction between 
the target and spectator nucleons, are also included to describe the MM £ j(7, 0) distribution 
of incoherent events 0, 3]. 

Furthermore, 0-meson events were singled out with a cut on the invariant mass of a 
K + K~ pair \M(K + K~) — M^| <0.01 GeV/c 2 , which is indicated as the region between two 
dashed lines in Fig. H^b) and Fig. Mb). For the jp — > <pp events, an additional cut on 
the missing mass |MM P (7, 0)| <0.03 GeV/c 2 is required. As seen in Figs. []Ja) and []Jd) 
and Figs. E]Ja) and[2^d), the contamination of A(1520) in the selected 0-meson events is 
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MM p (y,K + K") (GeV/c 2 ) M(pK") (GeV/c 2 ) 



FIG. 1. Spectra for the ^(7, K + K~)X reaction within the experimental acceptance: (a) Dalitz 
plot of K + K~p final state, (b) the invariant mass spectrum of K + K~ , (c) missing mass spectrum, 
and (d) the invariant mass spectrum of K~p for the selected <j) events. The dashed lines on the 
invariant mass of K + K~ in (a) and (b) label the region for the selection of <j) events. 

insignificant. 

The decay angular distributions of W{cos6), W((p), W($), and W(ip — $) in the helicity 
system for <ft events within the experimental acceptance are shown in Fig. [3] and Fig. HI 
The polar angle distribution W(cos9) behaves as ~(3/4) sin 2 9, indicating the dominance of 
s-channel-helicity-conserving processes. That the azimuthal angle ip of decay daughter K + 
aligns along that of the electric polarization vector $ of the incident photon signifies a larger 
contribution from natural-parity exchange processes. Such a correlation between the two 
angles, <p and $, is seen to be even stronger in the production from deuterium. Quantitative 
information will be given in Sec. IIVI 

Since not all particles in the final state are detected, there is a mixture of two compo- 
nents, incoherent and coherent, in the selected events from deuterium. To disentangle the 
individual contributions, extra effort is needed. By following the prescription specified in 
Refs. 

aa 

, the percentage of incoherent and coherent events, 1Z 1 and TZ 2 , are determined in 
the two separated missing mass regions, respectively, by a fit of the missing mass spectrum 
MM d (7, (p) with the MC-simulated distributions. The division of these two regions, MM div , 
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FIG. 2. Spectra for the d(7, K + K~)X reaction within the experimental acceptance: (a) Dalitz 
plot of K + K~p final state, (b) the invariant mass spectrum of K + K~ , (c) missing mass spectrum 
assuming rest deuteron target, and (d) the invariant mass spectrum of K~p assuming rest proton 
target for the selected <f> events. The dashed lines on the invariant mass of K + K~ in (a) and (b) 
label the region for the selection of <j> events. The MM^(7, <p) spectrum in (c) is fitted with the sum 
(solid line) of MC-simulated components of coherent (dotted line) and incoherent (dashed line) 
events. 

is chosen to be 1.89 GeV/c 2 . Event by event, relative weights composed of either 7Z 1 or 
1Z 2 are assigned to the likelihood of incoherent and coherent processes, depending on where 
the missing- mass of the event sits. In contrast to Eq. (12T]) . the likelihood function for each 
event is represented as the weighted sum of the individual likelihoods from the incoherent 
and coherent processes: 

-ln£(p inco ,p co ,iV) = 

- twm^i 1 *? + (i - n^;:^}} + n. (22) 

~^ wO nco ) v ; w(p co ) JJ 

Here TZi is either TZ l or TZ 2 according to the associated missing- mass value of the event; 
the spin-density matrix elements for the incoherent reaction p mco and coherent reaction p co , 
together with the expected number of events, N, are the parameters to be determined in 
the extended maximum likelihood fit. 
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FIG. 3. Decay angular distributions for the p(j, K + K )X reaction in the helicity system within 
the experimental acceptance: (a) TF(cos6>), (b) W(<p), (c) W(&), and (d) W(y> — <£). 

The extended maximum likelihood fit was performed in the framework of the MINUIT 
package. To locate the global minimum, the initial values of fit parameters were chosen to 
be the converged results from many test fits where a random point in the allowed range of 
the multidimensional parameter space was used as the start for the minimization. After a 
call of minimization by the method of MIGRAD, the statistic error estimation was done by 
a MINOS error analysis. The returned symmetric parabolic error was reported. 

Several procedures were checked to ensure the fit quality. At first, the fit status returned 
by MINUIT was required to have a normal convergence. The global correlation coefficient of 
each parameter should be greater than zero and less than 0.99, to avoid a wrong estimation of 
the statistic error for the case of uncorrelated or strongly correlated parameters. The typical 
value of the global correlation coefficient was distributed between 0.2 and 0.6. Second, the fit 
parameter for the number of events must be consistent with the input statistics, an advantage 
from the use of an extended maximum likelihood fit. Finally, reasonable agreement between 
the one-dimensional angular distributions of cos#, ip, $, and 93 — $ from the input events and 
the fitted values of p, was checked. The x 2 value per degree of freedom from a normalization 
fit is required to be in the range of 1.0-3.0. In general, a fit for a given kinematic bin 
associated with small statistics was likely to fail in this quality check. 
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FIG. 4. Decay angular distributions for the ^(7, K + K )X reaction in the helicity system within 
the experimental acceptance: (a) W(cos0), (b) W(tp), (c) W(3>), and (d) W(ip — <3?). 

The spin-density matrix elements do not depend on the beam normalization which is typ- 
ically the main systematic uncertainty of cross-section measurements at LEPS. We evaluated 
the systematic uncertainty in the following way. Ensembles of MC events were generated 
with specific sets of spin-density matrix elements being the same as those found in the 
real data. Those MC events were then filtered by the experimental detection efficiency. 
Subsequently, those events with statistics similar to the real data were analyzed for the de- 
termination of spin-density matrix elements p. The mean of the distribution of obtained p 
for many trials, the so-called pull distribution, was determined. The deviation of the mean 
from the generated value of p for MC events contributed to the estimation of the systematic 
uncertainty. 

Another source came from the background reactions [e.g., nonresonant K + K~ and 
A(1520) production]. This bias was estimated by comparing the results in the signal region 
[\M(K + K~) -M^| <0.01 GeV/c 2 ] and the sideband region [0.01 < \M(K + K~) -M^| <0.02 
GeV/c 2 ], by taking into account different signal-to-background ratios in these two regions . 
In general, the signal-to-background ratio in the defined 0-event region is good enough as 
shown in Figs. Ufa) and[T]^b) and Figs. |2^a) and|2]^b). Thus this bias was found to be small 
(of the order of 0.01—0.05 except at regions of larger \t\). Since the statistics in the sideband 
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region was not always enough for a reliable fit, we conservatively included this bias in the 
systematic error, rather than applying the corresponding correction to the results. 

For the 0-meson production from deuterium, there are additional sources of systematic 
error: the off-shell effect for the incoherent process in the MC simulation and the division 
point in the missing mass regions for disentanglement [2, 3]]. The choice of the division point 
MMdiv would affect the contents of coherent and incoherent events on the two separated 
regions. Four choices of MMdi v , 1.875, 1.88, 1.89 and 1.90 GeV/c 2 , were used for evaluating 
the systematic uncertainty. 



IV. RESULTS 

We measured the spin-density matrix elements of 0-meson photoproduction in the region 
of 1.57 < £ 7 < 2.37 GeV and \t\ < 0.2 GeV 2 /c 2 . We define t as t - f min or t - t d min , 
which is the squared four-momentum transfer t subtracted by its minimum value for the 
corresponding photon energy bin, under the assumption that a proton or deuteron is at rest. 
The binning sizes were 0.2 GeV for E 1 and 0.05 GeV 2 /c 2 for t. The real or imaginary parts 
of spin-density matrix elements pfj for the 7p — > (ftp, 'yd — > (ftpn, and 7c? — > (ftd reactions were 
determined by the method as described in Sees. [XT] and 



A. Production from free protons: "fp — > (j>p 

The spin-density matrix elements of the jp — > (ftp reaction in the helicity, Gottfried- 
Jackson, and Adair systems are shown in Figs. [U |6] and El respectively. Within the error 
there is no strong energy dependence in the measured region of 1.77 < < 2.37 GeV. Clear 
nonzero values for p\_ 1 and lmpf_ l are seen in all three systems, and p\_ 1 ~ 0.2. According 
to Eq. ([5]), a sizable 30% contribution of unnatural-parity exchange processes, other than 
the the natural-parity Pomeron exchange, is observed. 

In the helicity system, small but nonzero positive values for p[j and negative ones for 
Rep'j'o are observed. The p[j (R- e Pio) are °f positive (negative) values. The finiteness of Pq 
indicates the presence of amplitudes violating s-channel-helicity in the (ft production. 

In the Gottfried- Jackson system, Pq Q seems to be larger and Rep5 turns to be positive, 
compared to the results in the helicity system. The increase in the magnitude of p° suggests 
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FIG. 5. Spin-density matrix elements for the "yp —> <j>p reaction in the helicity system as a function 
of \t — t P i \ in various E~ regions. The vertical bars are for statistical error only. 
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a function of \t — t^-\ in various E~ regions. The vertical bars are for statistical error only. The 
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FIG. 7. Spin-density matrix elements for the jp —> 4>p reaction in the Adair system as a function 
of \t — t p min \ in various regions. The vertical bars are for statistical error only. 

a stronger violation of helicity conservation in the t channel. Specifically, those spin-density 
matrix elements reflecting helicity- nonconserving amplitudes such as Pq , p?_ 1; pj , Rep| , 
and Imp^ become clearly nonzero at large t. Figure |6] shows a comparison with the prediction 
from a model based on the dominance of the Donnachie-Landshoff Pomeron plus (tt, 77)- 



exchange channels |ll|] at E 1 = 2 GeV. This model correctly predicts the sign of each matrix 
element. However, the t dependence for p(] , and pJo i s n °t correctly predicted. Also 

the absolute scale of p\_ l and Impf_ x are overpredicted, which means that the contribution 
from the unnatural-parity exchange processes is not sufficient in this model. 

The p values in the Adair system are rather similar to those in the helicity system except 
for an opposite sign for Rep^o an d im Pio- 



B. Incoherent production from deuterons: ^yd — > 4>pn 

The spin-density matrix elements for the incoherent production from deuterons in three 
angular systems are shown in Figs. [HI M and [TU], respectively. The energy dependence is 
also insignificant. Because of statistical constraint, the measured range of t is limited to 
\t\ < 0.15. In general, the results are quite similar to those of the production from free 
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protons except that the absolute values of p\_\ and lmpf_ l are slightly larger (p\_ 1 ~ 0.25). 
This suggests that the contribution from unnatural-parity exchange processes is reduced in 
the production from neutrons. It could be interpreted as a destructive interference effect 
among the unnatural-parity (tt, ?y)-exchange processes in the 0-meson production from the 



neutron 
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111 



C. Coherent production from deuterons: jd —> <j)d 

Figures fTTl [12] and [13] display the spin-density matrix elements for the coherent production 
from deuterium in three angular systems. A distinct feature is a strong increase of the 
absolute values of p\_ x and Imp^_ 1; and p\_ 1 ~ 0.45, compared to what is observed in the 
reactions of jp — > <pp and 7c? — > <frpn. This again suggests that the contribution from the 
unnatural-parity exchange processes is significantly reduced. It could be understood as a 
result of the forbidden coupling of the isovector ir exchange with the isoscalar deuteron 
target. All the other components in p are similar to those in the production from nucleons. 
In the helicity system the smallness of pg (less than 10%) in the small- |i| region is the same 
as what was observed in the region of large \t\ by CLAS 
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FIG. 11. Spin-density matrix elements for the "yd —> <pd reaction in the helicity system as a function 
of Ej and \t — t^ml- The vertical bars are for statistical error only. 
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FIG. 12. Spin-density matrix elements for the — > <pd reaction in the Gottfried-Jackson system 
as a function of E~ and \t — tfL \. The vertical bars are for statistical error only. The solid lines 



are the theoretical predictions of p for the 7p 



reaction at EL = 2 GeV 
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FIG. 13. Spin-density matrix elements for the *yd — > <pd reaction in the Adair system as a function 
of E 1 and \t — t^ml- The vertical bars are for statistical error only. 



V. DISCUSSION 

The finiteness of the spin-density matrix elements p° , p?_ l5 Pq , Rep 10 , and Imp 2 is 
clearly observed in the Gottfried- Jackson system. This suggests the presence of helicity- 
nonconserving effects in the t-channel exchan ge p rocesses for 0-meson photoproduction, 



23j). Relatively, the helicity-nonconserving 



(e.g., two-gluon exchange as mentioned before 
effects become less in the helicity (s-channel) and Adair systems. 

Previously, we reported a nonzero value of p\_ Y (0.12 ± 0.03) for the jp — > (ftp reaction 
in the region of \t - t p mm \ < 0.2 GeV 2 /c 2 at £ 7 = 1.77 - 1.97 GeV [1], while p°_! became 
less (about 0.04 ± 0.02) at E 1 = 1.97 — 2.17 GeV. The measurement was done with a fit 
on the one-dimensional azimuthal angle distribution W(u?) in the Gottfried- Jackson system. 

n n 

In Refs. [13l . [14| , the authors argued that p°_i, which could only come from a spin-orbit 
interaction, must be close to zero near threshold and should monotonically increase with 
photon energy. However, the results were obtained in the region of \t — i^jj < 0.2 GeV 2 /c 2 , 
instead of at t = t p min . The effect of t dependence shown in Fig. |6] should be taken into 
account. 

Using the current data set from the hydrogen runs, we repeated the same analysis. The 
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results are shown in Fig. [TU Figures ITffi a) and ITffi b) display the VK(<p) distributions and 
the fitted p°_ 1 for the 7p — >■ 0p reaction in the regions of \t — t p min \ < 0.05 GeV 2 /c 2 and 
1^ ~~ tminl < 0-2 GeV 2 /c 2 at three bins. Obviously there is some t dependence. The 
current measurement at \t — t p min \ < 0.2 GeV 2 /c 2 is consistent with the previous results. 
At the smallest \t\ bin, p°-i i s about 0.05, which is consistent with the results shown in 
Fig. [6j Within the error, we could not draw any conclusion on the monotonic increase of 
Pi_ 1 with the photon energy but at least there is no large jump in the energy dependence. 
It is interesting to note that the pj-i value at E 1 = 1.97 — 2.17 GeV, where the peak of 
the bump structure in the 0-meson production cross section appears, seems systematically 
smaller than those in the neighboring energy bins. 



The measurements of the diffractive production at low energies are qualitatively con- 
sistent with the description of natural-parity Pomeron and unnatural-parity pseudoscalar 
exchanges in the t channel [12J. The latter term includes the ((f), T]) exchange and accounts 
for about 30% of the production cross section in the 7p — > <pp reaction. In the coherent 
production process, the isovector ir exchange is eliminated, and the natural-parity exchange 
becomes quite dominant. This observation reflects the fact that the ^-exchange channel is 
relatively small. A slight increase in the contribution for natural-parity exchange is seen in 
the incoherent production, being consistent with the theoretical prediction, if one takes into 



account the destructive effect of n- and //-exchanges in the 'jn —>• <pn reaction 
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A bump structure in the differential cross section at forward angles is observed at = 2 
GeV in both the production from free protons and the incoherent production from deuterons. 
The origin of this structure should be common for production from both protons and neu- 
trons. Thus it is unlikely that the interference effect from the (vr,-^) exchange is responsi- 
ble jl4j ]. Judging from the smallness of p except p\_ Y and Imp 2 _ x in the helicity frame, we 
see that the helicity-nonconserving amplitudes is rather limited for the s channel in the very 
forward direction. There is also no significant variation in the spin-density matrix elements 
across the bump region. Therefore we conclude that this bump structure is also unlikely to 



be caused by any nucleon resonance states containing large ss content 18|, |l9) unless there 
are some complicated interference effects. 
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FIG. 14. (a) Azimuthal angle distribution W(ip) for jp ->• (f>p at \t - f min \ < 0.05 GeV 2 /c 2 in 
the Gottfried-Jackson system. The solid curves are the fit to the data, (b) Same as (a) but at at 
|t-CJ<0.2 GeWc 2 . 



VI. SUMMARY 



In summary, we presented the measurement of decay angular distributions of mesons 
produced from protons and deuterons with linearly polarized photons at very forward angles 
at E 1 = 1.77 — 2.37 GeV. Nine spin-density matrix elements representing the angular distri- 
butions in the K + K~ decay mode are obtained simultaneously by the extended maximum 
likelihood fit in three different coordinate frames. There is no strong energy dependence 
observed in the measured energy region, although some t dependence is found. Small but 
finite helicity-nonconserving effects are seen in the t channel but they become less significant 
in the s channel. 

Based on the measurement of p\_ 1 and Imp 2 , unnatural-parity exchange processes are no 
longer negligible near threshold, even though the natural-parity Pomeron exchange process 
still dominates. The unnatural-parity processes are significantly reduced in the coherent 
production from deuterium. This fact can be interpreted to appear as the result of the 
suppression of isovector n exchange. The averaged contribution of the unnatural-parity 
processes in incoherent production decreases slightly and could be due to the destructive 
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interference of the (tt, ?7)-exchange processes in production from neutrons. 

A theoretical model based on the dominance of the Donnachie-Landshoff Pomeron plus 
(.^-exchange eb aI1 ne.s Q gl ves a reasonable prediction at S 7 = 2 GeV but obv.ons.y 
there is room for quantitative improvement. New and comprehensive information on spin 
observables available from this work should help to differentiate the theoretical models. 
Hopefully, this will lead to a better picture of the 0-meson photoproduction at low energies, 
especially to shed light on whether any exotic channel is necessary to account for the bump 
structure at E n ~ 2 GeV. 

A manifestation of an unnatural-parity exchange component in the photoproduction of 
4> mesons near threshold regions was already observed. New programs at LEPS and CLAS 
aim to perform measurements at = 2.5 — 3.5 GeV, and to study double polarization ob- 
servables with polarized proton and deuterium targets. Such experimental efforts shall bring 
further understanding of the appearance of bump structure in the 0-meson photoproduction 
near threshold. 
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